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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO.' 1326 

AIRFOIL IN SINUSOIDAL MOTION 
IN A PULSATING STREAM 
By J. Mayo Greenberg 


SUMMARY 


The forces and moments on a two-dimensional airfoil executing 
harmonic motions in a pulsating stream are derived on the basis of 
nonstationary incompressible potential -flow theory, with the 
inclusion of the effect of the continuous sheet of vortices shed 
from the trailing edge. An assumption as to the form of the wake 
is made with a certain degree of approximation. A comparison 
with previous work applicable only to the special case of a 
stationary airfoil is made by means of a numerical example and 
the excellent agreement obtained shows that the wake approximation 
is quite sufficient. The results obtained are expected to be 
useful, in considerations of forced vibrations and flutter of 
rotary -wing aircraft for which the lifting surfaces are in air 
streams of variable velocity. 


INTRODUCTION 


The problem of an airfoil in arbitrary motion in a pulsating 
Btream arises in connection with rotating blades in forward motion, 
for example, helicopter blades. Use of two-dimensional theory for 
this problem leads to deviations from reality with respect to the 
position of the wake, but, as most of the wake effects arise from 
that part near the airfoil, the error should not be serious. 

Restricted results, applicable to an airfoil at a fixed angle 
of attack in a pulsating stream, have been obtained by Isaacs 
(reference l) . In the present paper, under the same linearizing 
assumptions as are made in the derivation in reference 1 but with 
explicit consideration and simplification of the form of the wake 
extending from the rear of the airfoil, the methods of Theodorsen 
(reference 2) have been extended to obtain the forces on the airfoil 
not only at a fixed angle of attack but also in arbitrary motion. 

The result is in closed form as compared with the Fourier series 
result of reference 1. Inclusion of the effects of arbitrary motion 



2 


MCA TN No. 1326 


of the airfoil presents the possibility of application to forced 
vibrations and to flutter of helicopter blades. 


SYMBOLS 


a* 

a 

b 

a 

0 

h 

x* 

X 

t 

V 

1 } 

CD 

C(k) a F + iO 
Hx) 

P 

P 

cr 

P 


position of torsion axis of wing measured from center 

nondimensional position of torsion axis of wing 
measured from center 

half chord of wing 

fixed part of angle of attack measured clockwise 
from horizontal 

varying part of angle of attack 
vertical displacement, positive downward 
horizontal coordinate 
nondimensional coordinate 
time 

stream velocity 
frequency 

circular frequency (2jtu) 
reduced frequency 
C function (reference 2) 

Dirac delta function (reference 3 ) 
local, static pressure 
air density 

a number determining magnitude of stream pulsations 
perpendicular force 

pitching, moment about x » a, positive counter 
clockwise : 
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q> 

<p r 

u 

if 

af 

Subscripts : 

a 

P 

v 

h 

0 


noncirculatory velocity potential 

circulatory velocity potential 

strength of wake discontinuity 

phase angle with respect to stream pulsations 

element of vorticity 

fixed part of angle of attack 
varying' part of angle of attack 
stream velocity 
vertical displacement 

with 3 or h, maximum amplitude; with x, wake 
position; with v, steady part 


DERIVATION OF FORCES 


Method . - In a derivation parallel to that of reference 2, the 
forces due to the noncirculatory flow and to the effect of the 
wake are treated separately, the usual assumptions being made: 
incompressible potential flow, two-dimensional flat -plate airfoil, 
small oscillations, and plane wake extending from trailing edge to 
infinity . 

Noncirculatory force and moment .- Consider an airfoil of 
chord 2b at an angle of attack a + p with respect to a stream 
having velocity v directed to the right (fig. l) . Let the 
angle a be constant so that the only variation in angle of 
.attack is due to a variation in the angle p. If the airfoil is 
allowed to rotate about the point x’ = a’ with angular velocity 
P (the dot over a symbol denotes differentiation with respect to ‘ 
the time t) and to move downward with a velocity h, the velocity 
normal to the airfoil is 


e(x’) = -h - (a + P)v + p(a’ - x’) 
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The potential for such a normal velocity function is (see reference 2) 


cp » hh i/l - x 2 + v(a + P)b ^1 - x 2 + Pt 2 Qx - \/l - x* 


( 1 ) 


where 


x » 


Using the equation of motion for nonstationary flow 


~ * -v(£ + h£\ 

dt \p 2 


where 

p local static pressure 

p air density 

w local velocity 
and substituting 


W a V + ^ 

ox 


gives the pressure difference at the point x as 


( 2 ) 


Integration of this local pressure difference over the length of the 
airfoil gives as the total force 


P = -irpb 2 p. + v$ + v(a + 3) - bap 


(3) 
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The noncirculatory moment about the point ia a is 
M a * b 2 


J~ Ap(x - a) dx = -2pb 2 ~ a ) + 2P V ^/" <p dx 


which becomes 


M, 


a 


-«pb 2 j-v 2 (a + p) ~ bav(a + p) -vh + b 2 ^ + a 2 )^ 


(3a) 


Circulatory force and moment .- The velocity potential of an 
element of vorticity -At 1 ' at a position xg in the wake and its 
image AT distributed over the airfoil is (reference 2) 


Vxxr 


M +aT ,- 1 

"2jT tan 


yi-Tx 2 v^TT 




(*) 


The element -AP moves to the right' relative to the airfoil with a 
velocity v. Thus 


^xx 

__Q. 

at 


$2a 


Substituting this expression and — 

integrating the effect of the entire wake 
the force 


into equation (2) and 
on the airfoil yields 


P 


r 

-pvb J 


*0 




U dxr 


(5) 
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■where U d xq is the, ©lenient of vorticity AT at the point Xq. 
The circulatory moment which is obtained from 

M a a b 2 J*_ Ap(x - a) dx 


is, similarly, 




(5a) 


The Kutta condition requires that at the trailing edge of the 
plate the induced velocity he equal to Oj therefore, at the point 

X =5 1 



9 r = ”4 % ^ 

Introducing the potential from equation (l) results in 


sX°l/^T oaxo ’ i + T(a;fe> + ^' a ) (6) 


The type of stream velocity encountered by a helicopter blade 
section in forward flight is 


( iayb) 

A1 + ere / 


v = v Q \l + ae 


( 7 ) 
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where ~ is the period of pulsation of the stream. The assumption 

is made that a < 1 , "because if cr > 1 there would occur a reversal 
of flow which cannot "be treated "by existing methods. In usual practice 
a > 1 only for radial positions near the hub of a helicopter "blade. 

For harmonic motions of the airfoil 



P 0 e 


i(c»3t + 


i 

h « hp© 


(v + V 


> 




( 8 ) 


where - and are the periods of variation in angle of attack 

and in vertical displacement, respectively, and and •>}/“. are 
phase factors. ....... . . p . 

Equation (6) "becomes, after substitution of equations (f j and 

( 8 ), 
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It is sufficient to .satisfy the Kutta condition separately for each 
of the terns on the right-hand side of equation (9)j therefore, a 
sheet of discontinuity of the following form is set up: 


t i(® T t"Vo) 

w.r a s{”- =%)+ v ' 


+ Ujje 


l(cDpt-lj^Xg) iC^t-kjjXg) 


+ U h e 


+ U T +|3© 


|(cM<»(3)t-V(3 3J C)] 


( 10 ) 


where 5(oo - xq) is the Dirac delta function (reference 3), and 
where. 


v 

V " *0 


k .3S. 

P 3 *0 


*h 


*h* 

v 0 


k 


(civ + 


'v+3 


v 0 


Justification for assuming this form for the wake is given in the 
appendix. . , 

Combining equations (9) and (10) and equating terms of 
corresponding time variation yields 
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1 

2rt 


4 - *o> a *o ■- V- 


or r a « Sicv^a (the stationary circulation) 


1 r.jSSTTT 

ST 4 y^TlV tao-oToa 



'xq + 1" ' -’ilcpSQ 

TLe Ixq s* 


icOpPo' b \ 


-«) 


i% 

e 


i f°° .M) + 1 -i%xo 

274 ln^”i u ^ e ^0 53 


1 f i / x o + 1 „ ” iJ V+p z o 
2n'J 1 I/'xq - 1 U v+p® 


dx Q = av 0 P 0 e 




Introducing these relations and equation (10) Into equation ( 5 .) gives 
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The quotient of integrals occurring in this expression is the 
C function defined in reference 2 as 


A CO 


*o - 1!a b . 

e dx* 


C(fc) a 


/xx 2 - 1 


a 


STTT - iky o 


e dx r 
J 1 * *Q ” 1 C 


The force due to the wake is, therefore, 


f lea t 

P a -2jrpvb<jv0a + orT 0 osC(3c v .)e v 

L 


^PoHI " a ) + v o p o_ c ( k p)® 

Uo^u-%) 

+ io^hoCCkjjje 

+ crv 0 P 0 C(V(3^ e " -*J 




(11) 


Applying the same methods to equation (5a) as were applied to 
equation (5) gives the moment due to the wake as 

M a = -jrpvb 2 + v(a + p) + ShQ- - 

+ 2*pvb 2 (a + — ^ <jv 0 a + av 0 aC(k 7 )e ^ 

+ + V] 

, io> v t] 

+ hCfkh) + cv 0 pO(k v+ p)e ^ 


(11a) 
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Adding equations (3) and (11) gives for the total force on the 
airfoil 


P a -jtpb 2 [k + vP + v(a + P) - ahp] 

f io^t 

- 2rtpvb <v Q a + ovQaCCiCyJe 

+ [b(| - a)p + ^QpJc(kp) + hCffeh) 


crv 0 PC(l^ +p )e 


iOiyt 


( 12 ) 


Adding equations (3a) and (11a) gives for the total moment about 
the point x » a 


Me 


-itpb 2 jjrb^g - a^p - vba(a + P) - "bah + b 2 ^g + a 2 ^ pj 
+ 2itpvb 2 ^a + £>)|vo a + ovQaCi'k^e^ 


V 


(i - a )p 


p + v 0 P 


C(kp) 


+,hC(k h ) + ofVQpCf^pie 


iOJ v t 


(12a) 


Equations (12) and (12a) are applicable in any case of simple 
harmonic motions of the airfoil, regardless of the manner in which 
the motion arises. For instance, change in the angle p may be 

due either to a twisting of the blade or to the fact that the plane 

of the helicopter disk is inclined to the direction of motion. An 
obvious extension to equations (12) and (12a) to include the effect 
of an aileron can be made by applying the results of reference 2 
to the aileron potentials and forces in the same manner as the 
present application to the potentials and forces of the entire 
airfoil. 

If the variation in displacement or in angle of attack of the 

airfoil or in stream velocity is made up of a superposition of simple 
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harmonic motions, equations (12) and (12a) must he modified in a 
manner indicated by the derivation. 

Numerical example. - lor the purpose of checking the results 
of the pr e sent derivation and in particular the assumption as to 
the form of the wake, a numerical example is given, and a comparison 
is made with the example given in reference 1 for the case of the 
stationary airfoil. A more direct comparison cannot he simply 
obtained because the form of the wake is not given explicitly in 
reference 1. 

For P s h * 0, equation (12) becomes 

P « - Jtpb 2 va - ' 2irpvb Jj 0 a + v 0 craC(k v .)e ia V^J 
If the stream is undergoing sinusoidal pulsations of the form 
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where Lq is the steady lift due to Vq alone. From the foregoing 
equation, after the indicated real parts have teen extracted. 


P_ 

*0 



+ 



cos ciyfc + o(l + F) sin 03yt 


0% . ' 0% 

~2~ 008 + ~7f~ 


sin 2»yt 


The same values of the parameters k v and a as those chosen ty 
Isaacs as representative of typical helicopter practice are used; 
namely, ky = 0.0424 and a « 0.4, and values of the functions F 
and G are taken from a tabulation in reference 1. Then 


£-- = 1.074 - 0.0395 cos 0) y t + 0.768 


sin oyfc 


- 0.074 cos 2oyfc - O.OO96 sin 2o>yt 


The value obtained by Isaacs is 

•£- a 1.079 - 0.0376 cos tult't 0.770 Sin CDyt 
L 0 

- 0.079 cos 2ciyt - O.OO697 sin aoyt 


- 0.00061 cos 3ajyt - 0.0050 sin 3cuyt + , . . 


As can readily be seen, these two expressions for the lift are 
in very good agreement; therefore, the present assumption concerning 
the fom of the >?ake is quite reasonable. ' It is also to be noted 
that this agreement holds for relatively large values of a and 
relatively small values of the frequency. Even better results can 
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■be expected for higher values of the frequency of pulsation of the 
stream, as is shown in the appendix; 


CONCLUDING REMARKS 


The assumption of a sinusoidal wake form "behind an airfoil in 
a pulsating stream has been shown to be adequate for the derivation 
of the forces and moments. The forces and moments obtained are 
expected to be useful in considerations of forced vibrations and 
flutter of rotary-wing aircraft for which the lifting surfaces are 
in air streams of variable velocity. 


Langley Memorial Aeronautical Laboratory 

National Advi sory Committee for ' Aeronautic s 
Langley Field, Va,, September 24, 1946 
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APPENDIX 


FOBMULATION OF THE WAKE BEHIND AN AIRFOIL IN A 
PULSATING STREAM 


Use of Dirac delta f -unction . - In order to account for the 
stationary lift on the airfoil, that is, the lift due to the fixed 
part yq of the stream velocity and due to the fixed part a 
of the angle of attack it is necessary to assume that the vorticity 
which has "been shed at some distant past time has travelled to 
an infinite distance from the airfoil and is localized there. 

The delta function has the properties (reference 3 ) 



S(c 


x) dx *> 1 



f(x)S(c 


- x) dx at f(c) 


■where f(x) is an arbitrary function. The first of these properties, 
when applied to the first term on the right-hand side of equation (10), 
gives 



x) dx — r, 


a 


so that the condition that the total vorticity for the stationary 
case he preserved is satisfied. The second property has the 
effect of localizing this total vorticity at the point x = c 
where c is chosen to he infinite. 

S inusoidal part of wake . - Consider an airfoil moving hack and 
forth harmonically in a uniform stream having a velocity Vq. If 
ovq is the maximum velocity of the airfoil, then the velocity of 
the airfoil at any time t is given hy 

i(JDyt 


V A = crv 0 e 
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and the position of the trailing edge relative to the equilibrium 
position of the trailing edge is (fig- 2) 


^0 “ y a dt “ S| ® 7 Q " ( A1) 


The vortex element at the time t and at the distance 
from the equilibrium position of the trailing edge of the airfoil 
is the same as that vortex element which ms shed from the airfoil 
at the time t - T, where T is the time necessary for the vortex 
element to travel with velocity Vq from the trailing- edge position 
at time t - T to the position Therefore, 

T = % jg - Sgft - T)J ' (A2) 

Suppose that the vorticity at the trailing edge varies sinusoidally 
with the frequency u = Then 

Clt 

- AT 0 e tot 

But 


£T(£,t) « j&IX t -r] 

&0 


Therefore 


io>( t-T ) 

AT(g,t) = &r 0 e 




* AT 0 e 


(A3) 
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Referring now to a coordinate system the origin of which is at the 
center of the moving airfoil, that is, reverting to the case of the 
fixed airfoil in a pulsating stream, requires the transformation 


^ - i = g - g 0 (t) 


Equation (A3) "becomes 


AP(xQ,t) ** ATo© 


[xb-l4 0 ( t) -£ 0 ( t-T ) j 


(A^) 


From equation (Al) it can be seen that Cq 0 as . 

Therefore, as 00 


to 

AT(xo,t) ~ AT Q e 



AT1 i(cot-laco+k) 
AT q© 


(A5) 


where k = — . This relation shows that for large values of the 
v 0 

frequency of pulsations in the stream, the wake form approaches 
the sinusoidal form, which was used in the derivation of forces. 
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